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Abstract 

We examine the non-stationary evolution of a coupled qubit- 
transmission line-resonator system coupled to an external drive and 
the resonator environment. By solving the equation for a non-stationary 
resonator field, we determined the requirements for a single-shot non- 
destructive dispersive measurement of the phase qubit state. Reliable 
isolation of the qubit from the "electromagnetic environment" is nec- 
essary for a dispersive readout and can be achieved if the whole system 
interacts with the external fields only through the resonator that is 
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weakly coupled to the qubit. A set of inequalities involving resonator- 
qubit detuning and coupling parameter, the resonator leakage and 
the measurement time, together with the requirement of multi-photon 
outgoing flux is derived. It is shown, in particular, that a decrease of 
the measurement time requires an increase of the resonator leakage. 
This increase results in reducing the quality factor and decreasing the 
resolution of the resonator eigenfrequencies corresponding to different 
qubit states. The consistent character of the derived inequalities for 
two sets of experimental parameters is discussed. 

The obtained results will be useful for optimal design of experi- 
mental setups, parameters, and measurement protocols. 

1 Introduction 

The measurement of qubits is an important step in quantum information 
processing. For accurate qubit detection, the readout has to be faster than 
the qubit relaxation time. Moreover, to implement quantum error correction, 
the readout time must be less than the decoherence time. The relaxation pro- 
cesses can be of different origin (even due to the qubit- vacuum interaction). 
Also, the effects of measuring devices on the behavior of qubit states during 
the readout time are of great importance PQ-[3]. 

The simplest scheme for fast qubit readout of the phase qubit, based 
on the tunneling effect, (single-shot measurement) has been implemented 
and studied in |l]-[5]. According to this scheme, the measurement pulse 
adiabatically reduces the barrier between the potential wells, one of which 
forms the qubit states. As a result, the qubit in the upper state switches by 
tunneling into the neighboring well with probability close to one, while the 
qubit in the lower state remains unchanged. This kind of readout is limited 
by the strong current noise back action from the measurement device on the 
qubit. Moreover, this demolition measurement process destroys the qubit 
states which is unacceptable for most applications. 

The scheme of dispersive readout of qubit states has evident advantages. 
In what follows, we consider the case in which an external device probes the 
qubit state indirectly via a transmission line resonator (a two sided cavity) 
weakly coupled to the qubit. This scheme allows repeated measurements. 
The presence of the qubit-resonator coupling causes the resonator eigenfre- 
quency to be dependent on the qubit state. Hence, the number of photons 
in the resonator, being dependent on the probe field-resonator detuning, de- 
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pends on the qubit state. (This number is maximal for zero detuning.) For 
this reason, the field going out of the resonator contains information about 
the resonator state as well as about the qubit state. 

The amplitude of the transmitted field depends on the photon number in 
the resonator and on the resonator leakage. The increase of both quantities 
results in an increase of the transmitted signal that in general improves the 
signal-to-noise ratio. However, both of these factors decrease the qubit relax- 
ation time due to the qubit-resonator interaction. Moreover, resonators with 
high leakage (low quality factor) lose their ability to filter out those external 
noises which penetrate into the low-temperature region from the "external 
world". Also, for a relatively large number of photons in the resonator, the 
qubit-resonator interaction cannot be considered to be weak, and the advan- 
tages of the dispersive measurement, based on the perturbation approach, 
disappear. 

The optimal choice of both the readout strategy and the setup param- 
eters can be facilitated by means of theoretical description of the physical 
processes involved in the course of readout. In what follows, we use our recent 
approach [S] which makes it possible to obtain explicit expressions for the 
photon numbers and the measurement-induced relaxation rate of the qubit. 
We derive a set of inequalities which are required to provide a single-shot 
dispersive measurement, and we present examples of the qubit and resonator 
parameters. 

2 Hamiltonian and equations of motion 

As in our previous paper [5], we consider a flux-biased phase qubit embedded 
in a symmetric transmission line by a set of distributed inductances and 
capacitances. This is shown schematically in Fig. 1. The qubit is coupled 
to the center of the transmission line by the capacitance, C g , leading to 
their effective interaction strength, g. The qubit transition frequency, u q , is 
tunable by varying an external flux, $, shown schematically in Fig. 1. It is 
assumed that u q as well as the frequency of the external drive, u^, is close 
to one of the eigenfrequencies of the transmission line, u r . In this case, one 
can consider the transmission line as a resonator with a lumped inductance, 
L r , and a capacitance, C r . Then, the transmission line can be described 
in terms of a harmonic oscillator characterized by the single eigenfrequency, 
ui r — (-Zv r C r ) ^. 
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The total Hamiltonian is, 



H = -^hu q a z + hu r (a^a + ~) + X! hu n (blb n + + 

ihgcjy ^af — a^j + i% /„ ^6 n a* — &^a^ + i^/o ^ cc ^ — c ^ > (1) 

where cr^ is the Pauli operator, jb^c^a^njc) are the creation (annihila- 
tion) operators of the resonator, bath, and drive excitations, respectively. 
The bath is modeled by an infinite set of harmonic oscillators. (See Refs. 
[TO] and P].) 

The qubit-resonator interaction term (the fourth term in (pQ)) couples the 
variables of both subsystems. It is evident that only even harmonics of the 
transmission line are coupled with the qubit. We consider that the second 
harmonic has a frequency close to the qubit frequency. In the center of the 
transmission line, the voltage of the second harmonic is equal to pQ, 



V(t) = -^(a + a+), 

where L and c are, correspondingly, the length of the line and its capacitance 
per unit length. For a given term for the voltage, V(t), applied to both, C g 
and the qubit loop, we can easily express the interaction strength, g, as a 
function of other parameters of the system. This interaction strength, g, is 
given by the expression, 
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For calculational simplicity, we approximated the potential energy of the 
qubit by a parabolic dependence on the effective coordinate. More details 
can be found in |12j . 

As one can see from Eq. [H the interaction term is linear in the qubit and 
the resonator variables. Alternative interaction Hamiltonians (nonlinear in 
the resonator variables) describe the situation considered in [13] - [16]. The 
resonator-thermostat and resonator-drive interaction terms are similar to the 
resonator-qubit interaction. The analytical approach in Ref. [9] is based on 
a set of inequalities, 

Wi,Uj, Uij = CUi — LOj, (2) 
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Figure 1: A scheme of the qubit-transmission-line system. Cj is the capaci- 
tance shunting the Josephson junction, JJ. C g is the capacitance coupling the 
qubit and the transmission line. $ is the external flux threading the qubit 
loop. 

where the indices i,j indicate the corresponding frequencies, 0J q>r ,d- The 
relations (J2J) indicate that all frequencies are close to each other. 

Moreover, inequalities (T5]) were complemented by the requirement of a 
dispersive regime for our Hamiltonian. Usually, this is written in a form, 



« 1. 



(3) 
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Taking into account the fact that the interaction term depends linearly on a 
and a\ a more accurate criterion can be written: 
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where h r is the characteristic photon number in the resonator. 

Considering all of the operators in the Heisenberg representation, we have 
derived in [9] the equation of motion for the resonator field, a(t), 



d t + i(u r - x<?z + 



foe + J2 fnK - ig°i 



(5) 



where x = fl ,2 / w gr, and the operators in the right-hand side have the following 
time-dependence: 



c(t) = c(t )e-^-* o) , b n (t) = 6 n (t )e^ (t - to) , a+{t) = a + (* )e~^ ( <-H 
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The last term in the square brackets comes from the bath dynamics, and 
represents the linear damping of the cavity mode. The right-hand side of Eq. 
(j3J), in which operators evolve as if there is no interaction with the cavity, 
can be interpreted as the input mode. One can visualize that assuming that 
in the remote past (at to, t > to) the corresponding excitations were moving 
towards the cavity but were not still affected by it. 

In deriving Eq. (J5J), we used Eqs. (J2J) and (jlj. Also, o z was assumed as a 
slowly varying quantity. The consideration below follows the general scheme 
outlined in [10] and [TT] . 

Ignoring the initial conditions and the transient stage, the solution of Eq. 
(15]) can be written as: 



a(t) 



ifocjt) 
u dr + in/2 



ifnK(t) 



0O r . 



+ in/2 



+ - 

to 
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where u ir = oji — u r + \a z . 

Strictly speaking, there cannot be steady state of the system with the 
qubit in the excited state. Nevertheless, there are long-lived quasi-steady 
states with given values of a z = ±1 if the qubit relaxation time is much 
greater than that of the resonator (i.e. k^ 1 ). To understand at which pa- 
rameters and time intervals Eq. (jSJ) can adequately describe the state of the 
cavity field, an evolution equation for <j z is required. Within the dynamics 
governed by the Hamiltonian (JTJ, it follows from [9] that, 



dt 



a. 



2g 2 



qr 



2u qr sm(u qr t)e 
\foc\ 2 



+ K 1 



COS 



u qr t)e~ 



1 



1 



a z \n„-\ — 

2 V r 2 



2ojdr[ 1 — cos(uj qr t)e 2 ] + Ksm(u qr t)e 



(7) 



where n b r is the number of photons in the resonator, averaged over bath vari- 
ables. It is assumed that other relaxation mechanisms are not as important 
as those included in Eq. (JTJ). 

Considering only the initial stage of the excited state relaxation under 
resonant conditions (Cbdr = 0), we can omit the second term in the braces of 
(J7|). Then, Eq. (J7]) reduces to, 



9 f 2X , v re j 

777 1 &z H co%{u qr t)e 2 

at I oo qr 
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(8) 
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where the condition, \co qr \ » k/2, consistent with the dispersive regime, was 
assumed. By neglecting the right-side term, we obtain an equation describing 
the small-amplitude Rabi oscillations with frequency, u qr , in the vicinity of 
<j z — — 1. Its solution is given by: 

a z (t) « -1 + -^(l - cos(co qr t)e-^ (n b r + lj . (9) 

The oscillations decay for times greater than 2/k. Their influence on the 
measurement fidelity was analyzed in [9]. 

The total relaxation of the qubit can be described by taking into account 
the right-hand side of Eq. (JS}. It can be easily seen that the relaxation rate, 
7 r , is given by, 

7r = — 4 n r + 1 )- (10) 

Cdqr \ J 

As we see, the qubit relaxation occurs even in the case of empty resonator 
{n b r = 0) due to the resonator leakage, k. In the literature, this effect is 
known as a vacuum-induced relaxation. 

The condition for nondestructive measurement of the excited state can 
be expressed in the form of the inequality, 

7r (t-t )«l, (11) 
where, t — t , is the measurement time. 



3 Non-steady state of the resonator field 

Due to the measurement, the resonator state can vary with time. Following 
Ref. [32], we consider the dynamics of a resonator that was initially in 
the steady state corresponding to the qubit in the ground state. Only an 
insignificant number of photons is in the resonator because of the cavity- 
drive is detuned by —2%. It is assumed that at t = t , a strong 7r-pulse is 
applied directly to the qubit at its transition frequency. After this, the cavity- 
drive resonant conditions are satisfied (ua = u r + x), and the field amplitude, 
a(t), starts to grow. Using the approach of [9], we can analytically obtain an 
expression for the resonator field. With the initial (at t = to) value of a(t) 
given by steady-state term, Eq. fl6]), we can easily solve Eq. (jSJ), 

a(t) = e ( -^- K/2)( *-* o) a(t ) + f ^' e (-^-«/2)(t-t') x 
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foc(t') + ^2f n b n (t')-iga + {t') 



(12) 



Retaining only the first term in square brackets of Eq. (j!2p and integrating 
over t', we get: 



a(t) = f c(t) 



e -K(t-t )/2 i _ e -/t(t-t )/2-| 

-i2 X + k/2 + re/2 



(13) 



The first term in square brackets of Eq. ([TBI is due to the "memory" of 
the "steady state" . The second term describes the increase of the resonator 
field due to the qubit excitation. In the absence of the qubit transition, the 
resonator field would be: 

a(t) = /oC(t) , . (14) 

It can be seen from Eqs. ([TBI and (fT4"l) that the difference in the resonator 
fields corresponding to different qubit states vanishes for short measurement 
times, t — t <<<2/re. To notice the difference in outgoing signals, the 
conditions, 

«/2> (t-t )~\ (15) 

as well as, 

2 X > k/2 (16) 

must be satisfied. 

It follows from Eqs. (fl5l) and ([TBI) that 2\ > (t — to) -1 - This inequality 
can also be derived from the following general consideration. The spectral 
width of any (t — to)-lo n g pulse is not smaller than, (t — to) -1 - Therefore in 
the opposite case, 2\ « (t — to)" 1 , two harmonics with frequencies, u r ± x, 
being close to one another within this spectral interval, have approximately 
equal amplitudes. So, these harmonics cannot be easily resolved. 

There is an essential difference in the values given by Eqs. (fT3|) and (|T4|) . 
if t — to > 2/k. In this case, the value of outgoing field, b out (t), provides useful 
information about the qubit state. Using input-output theory [10]- [UJ and 
considering the resonator as a symmetric two-sided cavity, the transmitted 
field is determined by the resonator field as, 

bout(t) = -^|a(t). (17) 
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The average photon number leaving the resonator through the right side per 
unit time, is given by, 

(btut(t)bout(t)} = ^(a + (t)a(t)}. (18) 

Substituting a(t) from Eqs. f|T3|) or ffT4l) and considering c(t) as a classical 
variable, we can easily calculate the outgoing photon flux. 

The total number of photons, n^, leaving the resonator during the mea- 
surement time, t — t , is given by, 

n T = ^(t-t )n= £ f dt'{a + (t')a(t')). (19) 

It can be seen that by definition, n is the photon number in the resonator 
averaged over the measurement time, t — t . The required value of n can be 
achieved by varying the drive. For theoretical estimates, Eqs. (|T3|) . (|T4|) . and 
(TT9|) can be used. 

To get reliable information about the state of the qubit after the 7r-pulse 
(if it is in the ground or excited state) in the course of a single-shot measure- 
ment, many photons («t >> 1) must interact with the measurement device. 
It follows from Eq. ( fl9l) that this condition can easily be satisfied by increas- 
ing the measurement time and (or) by increasing the leakage parameter, k. 
Unfortunately, both ways have evident disadvantages in view of possible ap- 
plications of qubits: (a) quantum computers require fast measurements and 
(b) large values of k result in small quality factors, Q, for the resonators 
(Q = u r /n). The last circumstance can shorten the lifetime of the excited 
qubit state. 

The interplay of several factors should be taken into account for the op- 
timal choice of the measurement strategy. In the next Section, we consider 
this issue in more details. 

4 Conditions required for a nondemolition sing 
shot dispersive measurement 

We start from the inequality (flT|) . which provides a qubit to remain in the 
excited state (with high probability, close to unity) after the individual mea- 
surements. For the non-steady state, expression (TiTj) should be rewritten 
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as: 

JL-(t-t )n«l, (20) 

where n is defined by Eq. ( 1T9|) . In Eq. ( 1201 . we have neglected unity in 
comparison with the average photon number, n, in the resonator. Eq. ( 1201) 
can be rewritten in the form, 

f(t-to)n«^ (21) 

which is more convenient for physical interpretation. In particular, the left- 
hand side term is the total photon number which left the cavity during the 
time, t — to. Eq. ( 12T1) establishes an upper bound on this number. At 
the same time, this number should be sufficiently large. In this case, the 
measurement process acquires classical properties in spite of the fact that a 
quantum state is measured. Thus, we have, 

i«5<* -*>)»«$ = %■ < 22 > 

It seems that the second inequality can be easily satisfied by choosing 
a large qubit-resonator detuning, u qr , or a small qubit-resonator interaction 
parameter, g. However, the variations of these quantities are restricted by 
the requirement, 

21x1 = 2^- > k/2, (23) 

\LO qr \ 

which means that nonresonant eigenfrequency should be beyond the res- 
onator bandwidth. 

Inequalities (|22|) and fl2"31 should be complemented by two more which 
were discussed in the previous Sections. The first inequality concerns the 
validity of the dispersive approach, and follows directly from Eq. (j3J). In the 
case of non-steady state, it is given by, 



' X « 1. (24) 



\COqr\ 



The other inequality deals with the problem of the resolution of the two 
resonator eigenfrequencies during the measurement time, 

k/2> (t-t )-\ (25) 
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In summary, conditions (122]) - (1251) should be satisfied for nondestructive single- 
shot dispersive readout of the phase qubit. Are these conditions non-contradictory? 
Is it possible to satisfy all of them, at least in some particular cases? In the 
next Section we will illustrate such possibilities. 

5 Examples 

(i) We first choose the frequencies of the resonator and the qubit, and the 
coupling strength as: <jj r /2rr = 6.5 GHz, uj q /2-n = 8 GHz, and g = 0.05u qr . 
Then, the frequency shift is: xjlv: = 3.75 MHz. If the measurement time is 
a given parameter, for example, t — to = 4* 10~ 8 s, then the leakage value, k, 
can be taken as: 0.5 *10 8 s -1 . This quantity is sufficiently large to distinguish 
the contributions in b out arising from steady- and non-steady states of the 
resonator (see Eq. (TIB"]) ). At the same time, it is not large enough to violate 
the conditions fl2j) and ((23]). 

Finally, it is reasonable to assume = 10, hence, n = 10. The value of 
tit is sufficiently larger than the fluctuations of the outgoing photons (in the 

case of Poisson's statistics \J (Stit 2 ) = \fnr)- On the other hand, conditions 
( 122]) and ( l24j) are not violated in this case. 

The relaxation rate of the qubit due to the interaction with the resonator 
is: 7 r = 2.5 * 10 6 s -1 . Then, the probability of remaining in the excited state 
after the measurement, determined by 1 — j r (t — t ), is 0.9. 

The resonator quality factor is equal to Q = u r /n ~ 820. The relatively 
small value of the quality factor is due to the large leakage from the resonator. 
The power, coupled to the outside of the resonator, P = ^nhco r , determines 
the corresponding voltage, (V out ), via the relationship, P = (V out ) 2 /R, where 
R is the impedance at the output of the resonator. Usually R = 50 Q. Then, 
{V out } = 0.23fiV. 

(ii) For comparison, let us consider resonator and the qubit frequencies 
that are a factor of two larger than in the case (i): uo r /2ir = 13 GHz, uj q /2-K = 
16 GHz. 

As before, g = 0.05u qr . Then x/2vr = 7.5 MHz. If the photon numbers 
and measurement time remain unchanged (n^ = 10, n = 10, t — to = 4 * 
10 _8 s), then k = 0.5 * 10 s s _1 . Because of the increase of u r , the quality 
factor as well as the ratio x/ K doubles. A similar tendency concerns the 
increase of the measurement time. For greater t — to, the value of k can be 
smaller resulting in a larger quality factor. Also, in this case, the output 
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voltage, (V out ), increases by a factor of \f2. 

The increase of the frequencies makes it easier to choose the desirable 
parameters required for a single-shot dispersive measurement. For example, 
it becomes possible to decrease the relaxation rate, j r , improving in this way 
the measurement fidelity. To show this, let us assume in the case (ii) that 
g/co qr = 0.05/\/2. Then, the relaxation time of the excited state doubles, and 
the measurement fidelity grows from 90% to 95%. In spite of the decrease of 
x/ k toward the value calculated for the case (i), the resonator lines remain 
well- resolved. It is important to emphasize, that a similar variation of g/uj qr 
in the case (i) will result in a very small x being insufficient for resolving the 
lines. 

6 Conclusion 

At first sight, the statement that a single-shot dispersive measurement of 
the phase qubit can also be nondemolition measurement appears unrealistic. 
A simple analysis of inequalities (|22l)- (l25l) shows that it is rather difficult 
to make a nondemolition measurement. Nevertheless, we have shown, for 
particular cases, that this measurement is quite realizable, although only for a 
restricted range of the parameters. For performing this kind of measurement, 
it is important to take into account the fact that shortening the measurement 
time requires a growth of the leakage, k, thus decreasing the resolution of the 
resonator lines and decreasing the quality factor. Moreover, a very strong 
drive can violate the conditions required for the dispersive measurement and 
it can shorten the lifetime of the qubit in the excited state. On the other 
hand, using resonators and qubits with high frequencies enables the choice 
of suitable parameters. 

In this paper, we have only touched the problems of the qubit-resonator- 
drive noises that are very important for quantum measurements. The mea- 
surement induced noise arises from the second and third terms in square 
brackets of Eq. ffl2j) . They can be accounted for straightforwardly. Never- 
theless, the complete solution of this problem is possible only by accounting 
for the measuring scheme and the intrinsic noises of the measuring device, 
and by taking into consideration the noise produced by the amplifiers. All 
these issues require further consideration. 
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